Field theories and exact stochastic equations for interacting particle systems 
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We present a new approach to the dynamics of interacting particles with reaction and diffusion. 
Starting from the underlying discrete stochastic jump process we derive a general field theory de- 
scribing the dynamics of the density field, which we relate to an exact stochastic equation on the 
density field. We show how our field theory maps onto the original Doi-Peliti formalism, allow- 
ing us to clarify further the issue of the 'imaginary' Langevin noise that appears in the context of 
reaction/diffusion processes. Our procedure applies to a wide class of problems and is related to 
large deviation functional techniques developed recently to describe fluctuations of non-equilibrium 
systems in the hydrodynamic limit. 
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Many problems of current interest in statistical physics 
involve strongly interacting particles that exhibit non 
trivial collective phenomena. One example is that of su- 
percooled liquids, where the dynamics slows down dra- 
matically as the glass transition is approached, due to 
the increasingly collective nature of the dynamics Q. 
Other examples are given by systems of diffusing parti- 
cles that branch and/or annihilate and that are subjected 
to an external drive; depending on the relative strength 
of these effects a variety of non-equilibrium transitions 
and anomalous scaling behaviour can appear 0. Devel- 
oping theoretical techniques for such difficult problems 
is of great importance, in view of the diversity of situa- 
tions in which they can apply, and of much broader scope 
than usual equilibrium statistical mechanics. A natural 
framework in this context is that of field theory which 
combined with perturbative renormalization group tech- 
niques have been applied successfully to a wide range 
of non-equilibrium phase transitions Furthermore, 
even in non perturbative regimes, field theory can be 
very useful because it allows one to articulate different 
types of approximations, such as the Exact Renormali- 
sation Group approach [J| and Mode Coupling Theory 

1^ . A field theoretical formulation of interacting par- 
ticle systems which has become somewhat standard in 
recent years is provided by the Doi-Peliti formalism (dp) 
. Starting from a second quantization representation of 
the Master equation, one obtains, after a rather elaborate 
coherent state representation, a field theory in terms of 
two fields (j) and (j) (see below) which has been the starting 
point of a very large number of studies However, be- 
sides its intrinsic complexity, the formalism appears to be 
fraught with difficulties. For example, the action of the 
field theory corresponds to a reasonable looking Langevin 
equation for the density of particles, except that the noise 
is often complex or even pure imaginary. This suggests 
that the field (j), despite its superficial resemblance with 
the density, in fact lacks direct physical interpretation, 
see e.g. [g. 

The aim of this Letter is to formulate a two-field the- 
ory for interacting particles in a transparent way, starting 



from a natural representation of the microscopic stochas- 
tic dynamics of the system in discrete space. Our for- 
malism focuses on the physical density field p and, as a 
consequence, is related directly to the stochastic equa- 
tions governing the evolution of p. The original dp for- 
malism is recovered performing a canonical Cole-Hopf 
transformation. Our method is straightforward, free of 
the ambiguities related to the imaginary noise and ap- 
plies to many different situations (diffusion in an ex- 
ternal force field, pair- wise interacting particles, branch- 
ing/annihilation processes, hard-core particles etc.). Fur- 
thermore it is related to recent techniques developed in 
mathematical physics (see below) 0. As an interesting 
physical side product, our field theory in the case of liq- 
uids is found to be identical to the one obtained from 
a Martin-Siggia-Rose-De Dominicis-Janssen (msrdj) 
representation of the stochastic equation on the density 
field derived by Dean for Langevin particles ^3 ■ This is 
important since Dean's derivation contains several sub- 
tleties. In particular, some have raised concerns about 
a possible hidden coarse-graining procedure that would 
explain why the ideal gas entropy appears, quite unex- 
pectedly, in the Langevin equation. We end the paper 
with various technical comments, and possible applica- 
tions and extensions of our formalism. 

Let us start from the simplest situation - two sites la- 
beled 1 and 2 between which particles hop back and forth 
with a Poisson rate W12 and W21. The (integer) number 
of particles on the two sites are ni and 71.2. The variation 
of Ui between t and t + dt will be noted dJi as is standard 
for Poisson jump processes it is not a small quantity 
since it is equal to or ±1, but the probability for it to 
be non zero is of order dt. Of course, dJi and dJ2 are 
strongly correlated since a particle leaving site 1 lands on 
site 2 and vice-versa. More precisely, dJi = —dJ2 = +1 
with probability n2W2idt, dJi = — dJ2 = — 1 with prob- 
ability niWudt, and dJi = dJ2 = otherwise. We now 
introduce, a la msrdj, the generating function for the 
histories of the system: Z({n, n}) = {Y[^ exp[ni(t)(dJi — 
dni) + iT.2{t){dJ2 — dn2)]), where the averaging (...) is 
over the realizations of the Poisson jump processes. From 
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the above rules, it is easy to find the result in the limit dt 0: 



Z({n,n}) = exp / dt [-hidtni ~ h2dtn2 + niWi2ie"'-"' - 1) + n2VF2i(e"i-"^ - 1)] . 

I 



(1) 



Obviously, one could add different processes, such as for 
example on-site annihilation, where dJj = — 1 with prob- 
ability riifidt^ branching, where dJi — +1 with probabil- 
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ity riiiydt, two-body annihilation where dJ, = —2 with 
probability 7ii(?ii — l)Xdt, etc. With iV sites on a lattice, 
the total MSRDJ 'action' reads S = \nZ: 



S{{n, h}) = f dth + E ".W-,,(e"-"' - 1) + E {/^(e""' - 1) + He"' - 1) + A(n, - l)(e-2"' - 1)}], 



where (ij) means that the sum is over all couples 
In addition, a factorized initial condition with distribu- 
tion p(rii{0)) can be included by adding a contribution 
EJ-n.(0)n,(0) + lng(n,(0))], with gix) = E,P('Z)e^"- 
As it is customary in procedures a la MSRDJ the aver- 
age value over the stochastic dynamics of any observable 
O, a generic function of {ni{t)}, equals the average of O 
over the field theory characterized by the action |2Il and 
the fields rii^fii which should be treated as continuous. 
This ends the derivation of the field theory that, as antic- 
ipated, turns out to be much more straightforward than 
the DP one. 
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We now want to consider the continuum limit of 
©. In this aim we write W^j = Wocxp{[Ui - Uj]/2T) 
where Ui is an on-site potential (possibly time depen- 
dent), which varies on scales much larger than the lat- 
tice spacing a and T is the temperature. Restrict- 
ing to nearest-neighbour hopping, we therefore write: 
Wij = Wo (1 — o-Bij ■ VU {x, t) /2T) , where x is the position 
in space of site i and is the unit vector pointing from 
i to j. Defining the local density field p{x, t) = Ui/a'^ and 
p{x, t) ~ hi and expanding to second-order in gradients, 
we finally obtain (in the a — > limit, and with A = for 
the time being - see below): 



Si{p. p}) = j dtdx [-p{x, t) dtp{3:, t) + p{x, t) {pie-P - 1) + i^ief - 1))] 



-7 / dtdxVp{x,t) ■ p{x,t)VUix,t) +-fT J dtdx [-V p{x,t) ■ V p{x,t) + p{x,t){V p{x,t))'^] 
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(3) 



where 7 = Woa^/2T is the mobility of the particles. Note 
that the above derivation is easily generalized to many 
other processes and is independent of the particular form 
of Wij as long as detailed balance is verified (although 
one could also consider non-potential force fields as well). 

Before mapping this action onto the more standard dp 
form, we want to specialize to the case where p = u = 0, 
and U{x,t) comes from a two-body interaction between 
the particles, i.e: U{x,t) = J dyp{y,t)V{x — y). In this 
case, one finds that the above action is identical to the 
one obtained by using the standard msrdj representation 
for the following Langcvin equation, derived by Dean for 



interacting Brownian particles: 



dtpix,t) = 



7p(f , t)V 



dF[p] 
dp 



V P(^> t)ff{x, t) 



(4) 



where rf is a Gaussian white noise with correlations 
{r]c,ix,t)r]0{x',t')) = 2-fTSaf3d{x-x')6(t-t') and F[p] is 
the effective free-energy which has the naively expected 
shape F[p] J dxp{x)\np{x) + 1/2 J dxdyp{x)V{x - 
y)p{y) [Note that Dean set the mobility 7 of the particles 
to unity] . This is quite remarkable, since p{x, t) in Dean's 
equation is the exact continuum microscopic density of 
the system, p{x,t) = ^i5{x — ri{t)), before any coarse- 
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graining [ri{t) denotes the position of particle i at time 
t]; it is therefore not at all trivial that the effective free- 
energy should have a mean-field form. The fact that our 
derivation, which starts on a lattice and never uses Ito 
calculus, leads to the same action confirms the validity 
of Dean's analysis. 

Let us now consider the following Cole-Hopf change 
of variables, from the above fields p, p to new fields (f), (p 
defined as: 

p{x,t) ^ (j>{x,t)(j>{x,t^); /5(f, i) = ln(0(f,t)). (5) 

The Jacobian of this transformation is 1 and thus the 
measure is preserved. The reason for ^ t — e will 
be discussed in the following. Note that this change of 
variable could also be defined on the discrete lattice and 
it was already considered in the literature 0, |^ , 
but the precise connection that we establish here seems 
not to have been noted before. It is straightforward to 
show that the above action S{{p,p}) transforms into: 

Si{(j)J})= J dtdx {^[-dtcj) + iTV^^ + jV{VU(j))] 

(6) 

Setting U = in the above equation, we recover ex- 
actly the DP action for the problem of diffusing, branch- 
ing and annihilating particle, usually derived in a rather 
thorny way from a second-quantization representation 
of the Master equation (up to some boundary terms, 
which for clarity will be discussed later). Note that for 
/i = z/ = 0, the above action can be seen as the two-field 
representation of the propagator of the non-Hermitian 
Fokker-Planck operator for particles diffusing in a po- 
tential field In addition, in the case of two-body 
interaction between the particles, the V[VC/^] term be- 
comes V[/ dy(l){y,t)(j){y,t)\/V{y — x)(j){x,t)]. Let us fi- 
nally remark that the mapping from the action (j^J to its 
DP counterpart works also on the lattice, the continuum 
limit does not play any important role from this perspec- 
tive. 

The above action looks very close to a msrdj rep- 
resentation of the naive Langevin equation describing 
the problem if one interprets <^ as a density and shifts 
— > -I- 1 in order to have {(j)) = as in usual msrdj 
(see e.g. 0). However this interpretation is problem- 
atic since, as mentioned in the introduction, the noise 
term is unphysical 0. For example, in the case of dif- 
fusing particles with pair-wise interaction (fj, = v ~ 0) 
the "noise" term has a correlator {r]{x,t)i]{x' ,t')) = 
'2'^x[4>{x){'VxV){x — x')(j){x')]. For a uniform density (j) 
and V{x) = Vb exp(— x^), the variance of the noise is a 
negative definite operator meaning that the noise has an 
imaginary part! On the other hand, we clearly see that 
this difficulty disappears when one uses the fields p, p, 



which encode a very well defined underlying Langevin 
equation, albeit with non Gaussian, Poisson jump terms 
corresponding to particle creation and annihilation. 

Let us now focus on some technical but important 
subtleties of the transformation (jSJ. Consider pair- 
annihilation [A + A ^ 0) with rate A on a single site. 
The corresponding contribution in the MSRDJ and dp ac- 
tions are \p{p — 1) (e"^'' — l) and A(/)^(l — (f)'^), which do 
not transform exactly into one another under Q. The 
underlying reason is that the two field theories corre- 
spond to different time discretizations. This makes a dif- 
ference when the action contains non linear terms evalu- 
ated at the same position in time and space, because the 
response function {4>{x, t)(f)(x, t')) is discontinuous when 
t ^ t' . This can be traced back to the fact that the 
fields (j) and (j) are the coherent state representations of 
creation/annihilation operators of the theory, a and a"*", 
which satisfy [a, a+] = 1. In the dp formalism the ac- 
tion is obtained after normal ordering Q and as a conse- 
quence, (j){x,t)4>{x,t)4>{x,t)4>{xTt) is in fact the continu- 
ous time limit of (l>{x,t)(j){x,t — e)(j){x, t — 2e)(j){x, t — 3e). 
Instead, the field theory obtained from through the 
transformation ((SJ has no normal ordering and there- 
fore the term Xp^ (af, t) is the continuous time limit of 
\4>{x, t)4>(x, t — t)4>(x, t — 2e)(j){x, t — 3e). Thus, to trans- 
form our field theory into DP, one has to use the trans- 
formation 10) and in addition to take care of the time 
discretization which amounts to do normal ordering, and 
to recover the —1 missing in the above example. An 
alternative way to make the exact connection between 
the two field theory is through the operator formalism. 
After having expressed the Master equation in terms of 
the operators a, a"*" one does a canonical transformation 
a = e^'' p, a+ = ef such that [p, P'^ = 1- The op- 
erators p,p^, originally introduced in |9|, lead to a dif- 
ferent representation of the Master equation. This then 
leads exactly, using again a coherent state representation, 
to the field theory that we derived above directly from 
underlying the stochastic process, including the correct 
boundary terms. 

Another case where our strategy applies is interact- 
ing particles evolving under Newtonian dynamics, for 
which Doi also derived a field theory Our previ- 

ous starting point consisted in deriving exact stochas- 
tic equations for the density field and then applying the 
MSRDJ procedure to get the field theory. For Hamiltonian 
dynamics, one can derive an exact deterministic equa- 
tion on the density in position r and momentum p space, 
p{x,p;t) = J2iSi'^-f'i)^i'P~'Pi)- 

dtp{x,p;t) — -p- d.sp{x,p;t) 

r (7) 
-I- / dx'dp'p{x\p';t)d3V{x - x') ■ dj;p{x,p;t). 

In this case only the initial conditions are stochastic. 
The MSRDJ field theory corresponding to this equation 
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maps exactly onto Doi's using the transformation, akin 
to 10), (j){x,p;t) = e~P^^'^'-*-^ p{x,p;t) and ${x,p;t) = 
g,p{x,p;t) (^jjQ^g ^}jat p{x,p;t) is the MSRDJ field conjugate 
to p{x,p;t)). 

Finally, our formalism in the continuum limit bridges 
the gap between the purely microscopic Fokker-Planck 
evolution operator and the hydrodynamics description 
studied recently in the mathematical physics literature. 
A method to handle both the hydrodynamic limit and 
large and rare fluctuations around it, has been developed 
in jg] ■ The starting point of this work is very close to the 
generating functional we used. However, in order to focus 
on hydrodynamics length-scales and timcscalcs, as done 
in one has to consider a conjugated field fii = p{x/L) 
that is constrained to varying only on length-scales of the 
order of the linear system size L. In this case the "action" 
or functional (0) becomes a function of the hydrodynamic 
density field that represents the average density inside 
very large boxes, and can be related to the rate func- 
tional introduced and studied in Q . On the other hand, 
the stochastic equations corresponding to our continuum 
limit arc valid on a scale £ much larger than the lattice 
spacing a but much smaller than the system size L. As 
a consequence these allow one to tackle, with field theo- 
retical techniques, dynamic phase transitions where the 
physically relevant length-scales are much larger than a 
but not necessarily much larger than the (diverging) cor- 
relation length ^. The hydrodynamic limit of '8| instead 
corresponds to length-scales much larger than ^. 

In conclusion our procedure allows one to derive rather 
straightforwardly a field theory different from, but dual 
of, the Doi-Peliti formalism. This could be useful in 
cases where this by now standard framework docs not 
work, see e.g. It certainly avoids, contrarily to dp, 

very cumbersome computations due to normal ordering 
in cases in which the rates are complicated functions of 
the local density. Another advantage of our approach 
is that it is, almost by construction, directly related to 
stochastic equations on the density field. The represen- 
tation in terms of a stochastic equation, especially after 
having taken the continuum limit, is particularly appeal- 
ing. It can be helpful for numerical investigations [l^ 
since it might be more efficient to integrate numerically 
than to simulate the original lattice model. Furthermore, 
stochastic equations are very useful to encode and study 
universality classes as it has been understood in the case 
of critical slowing down |l7l | and non-equilibrium phe- 
nomena such as surface growth Some applications 
of the results of this letter are underway, for example the 
study of the condensation phase transition in the zero 
range process jl9| . From a more general perspective our 
results make clear that exact stochastic equations can 
be always obtained, thereby avoiding phenomenological 
guesses which, especially for off-equilibrium cases, can be 



very tricky |20l |. 
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